Latent heat in the chiral phase transition 
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The chiral phase transition at finite temperature and density is discussed in the framework of 
the QCD-like gauge field theory. The thermodynamical potential is investigated using a variational 
approach. Latent heat generated in the first-order phase transition is calculated. It is found that 
the latent heat is enhanced near the tricritical point and is more than several hundred MeV per 
quark. 
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I. INTRODUCTION 

A many-particle system in general possesses two or 
more different forms of existence called phases. Conver- 
sions between these phases are called phase transitions. 
There are mainly two kinds of phase transitions. One is 
the first-order phase transition, the other is the second- 
order one. In the first-order phase transition, as in the 
conversion from vapor to liquid, the entropy of the system 
changes discontinuously and a certain amount of heat 
may be generated as a result. This kind of heat has been 
called latent heat since named by Joseph Black (1728- 
1799). On the other hand, in the second-order phase 
transition, the entropy changes continuously and the la- 
tent heat is not given off. 

The behavior of the phase transition in quark mat- 
ter is of fundamental interest pj. It is believed that 
at sufficiently high temperature and/or density the sys- 
tem lies in the chirally symmetric phase. Quark-gluon 
plasma (QGP) is one of the most important possibilities 
in this phase. Recently there has been great interest in 
generating QGP in the experiments. It is expected that 
QGP may be realized in high energy heavy-ion collisions 
at CERN Large Hadron Collider(LHC) and at the BNL 
Relativistic Heavy Ion Collider(RHIC). When the tem- 
perature of QGP falls and reaches the critical temper- 
ature (T c ), the chiral symmetry in QCD is dynamically 
broken and we have a phase transition from the chirally 
symmetric phase into the chiral symmetry broken phase. 

The chiral phase transition at finite temperature 
and /or finite density has been studied by many authors 
• According to them, it was shown that the transi- 
tion is of second-order at low density and of first-order at 
high density. The critical density between them is called 
a tricritical point • This nature of the phase tran- 

sition was shown by the NJL model Qpl and then con- 
firmed also by the QCD-like theory pj|-|lC|. If the phase 
transition is of first-order, a latent heat must be gener- 
ated across the critical temperature as mentioned above. 
Our main goal here is to evaluate the latent heat in the 
chiral phase transition at finite temperature and at finite 
chemical potential with the use of the QCD-like theory. 



So far much efforts have been made to study about sig- 
nals of generation of QGP 01 • The generated latent heat 
may give a signal of QGP. For example, the shock wave 
can be expected to be produced by the generated latent 
heat and may be available as a signal. The shock wave is 
supposed to be given off at the phase transition between 
the two phases. This process is due to the acceleration of 
particles by the generated latent heat. Therefore, knowl- 
edge about the latent heat is an important prerequisite 
for analysis of the data in such a high energy collision 
experiment that is expected to generate QGP. 

The outline of the paper is as follows. In the next sec- 
tion, we derive the partition function of the quark mat- 
ter based on the QCD-like gauge field theory in Landau 
gauge. In sect.V, we perform a mean field approximation 
and derive the Schwinger-Dyson equation for the effec- 
tive mass. Then (sect.) using a variational procedure, 
we solve numerically the Schwinger-Dyson equation with 
a trial mass function and determine the critical point. In 
sect.X, we compute the latent heat numerically. Sect.Y 
is devoted to summary and discussions. 



II. PARTITION FUNCTION 

We derive here the effective partition function of the 
quark matter in the framework of the QCD-like gauge 
field theory. In this paper, we take N c = 3 as the number 
of color and Nf = 3 as that of flavor. Let us start with 
the exact QCD partition function in the path-integral 
formalism |15lllq : 



Z(T,p) = Tr e~ m 
= N DVDVDAexp 







dr Id ie£qcd 
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, (1) 



where = T , t = it and H stands for the Hamiltonian 
of the system. When the chiral symmetry is exact, the 
QCD Lagrangian Cqcd is given by 

£qcd = C q + Cf, (2) 

where Lagrangians of the quarks and gluons are given by 

C q = W( l7 ^ + M7 °)*-.g*7^*^, (3) 
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nonlinear terms, 



(4) 



respectively. Here fi, g, X a and a denote chemical po- 
tential, the coupling constant, color SU(3) matrices and 
gauge parameter, respectively. The function D~l = 
g^iud^d^ — (1 — a~ 1 )d fl d v is the inverse tree-level gluon 
propagater. 

Next, we neglect the nonlinear terms and the coupling 
constant is replaced by the run ning coupling one, which 
is called QCD-like theory [13 -[HT Then we can inte- 
grate out the gluon held A. Moreover to investigate the 
antiquark-quark pairing correlation, we make use of the 
Fierz rearrangement to Cqcd- If the most attractive 
term is left, the Lagrangian is expressed by 



-QCD 

,9' 2 



d%{^{x)^{y)}D F {x - j/){*(y)¥(aO}, (5) 



where g' = (3 + a)g 2 /27 and the inverse quark propaga- 
tor G _1 (x) = {i^dfj, + ^7°) is introduced. We hereafter 
fix the gauge parameter as a = (Landau gauge), which 
is the usual choice in the QCD-like theory. 

Let us now introduce an auxiliary bilocal scalar field 
tp(x, y) representing a wave function of an antiquark- 
quark pair. In order to use the Stratonovich-Haverd 
transformation, we take the following identity: 



1 = C I Dip*D<pexp 



d A xd A y{ip{x,y)-^{x)^{y)} 



xg l2 D F (x-y){tp*(x,y) - V(y)V(x)} 



(6) 



where C is a normalization constant. Substituting this 
identity into Eq.(l), the partition function can be rewrit- 
ten as 



Z = N' J D^D^Dtp*Dipexp[-S'} 
where the corresponding action is given by 



(7) 



S' 



d xd y 



— \<p(x,y)\ 2 D F (x - y) 



-*(x){S(x - y)G-\y) + y —ip*{x,y)D F {x - y) 



—2 D f{v - x)cp(y, x)}V(y) 



(8) 



Here we assume that the auxiliary field is dependent 
only on the relative coordinates: <p{x, y) = ip(x-y). This 
means that the center of mass of the quark-antiquark pair 
is stationary, because we concentrate on the ground state 
of the system. Let us define the pairing field A(x — y): 



A ( x -y) = g' d f(x - y)<p(y, x). 



(9) 



Using Fourier transformation, the modified gluon propa- 
gator D F (x — y) is expressed as 



Dpix-y)^?- 1 ^ J 



(2tt) 3 -p 2 



(io) 



where p = (iu) n ,p), x = (— ir,x) and Yl! n denotes the 
boson Matsubara frequency sum over uj n — 2nirf3~ 1 with 
an integer n. In the same way, the auxiliary field is ex- 
panded as 

<p(x -2/)=r 1 E / ^<p(p)e~^y~*\ (11) 

where ^ denotes the fermion Matsubara frequency sum 
over uj n = 2(n + l)7r/3 _1 with an integer n. Then we can 
rewrite the pairing field as 

A(x J ^A(p)e-*(*-), (12) 

where the Fourier component in the right hand side is 
defined by 



d 3 q g' 2 



Ml)- (13) 



) 3 -(p- q) [ 



Using the A(p), we can transform the partition func- 
tion. After integrating out the quark fields $ and \&, we 
obtain 

Z =N'fD<p*D<pexp[-S eB ], (14) 
where the effective action S s is given as 

S cS = -2^]T J -0^\og{A 2 {p)-p 2 ) 



>E/|f,^)A(,), 



(15) 



where p = (iuj n + fj,, p) and £ = N c x Nj = 3 x 3 = 9. 
Thus our effective action is represented by only the A(p). 
It is interpreted as an effective mass for the quark so that 
we call it the effective mass function. The effective action 
obtained above will play a fundamental role in the next 
section. 



III. MEAN FIELD APPROXIMATION 

In this section, we will approximate the partition func- 
tion obtained above. To this end, we make use of the 
stationary (WKB) approximation for the path integral 
[TH IT^ . If the extremum of the integral is realized by 
ip , it must satisfy the following stationary condition: 



SS c s 



6ip(p) 



= 0. 



(16) 



3 



The is a mean auxiliary field and has physical mean- 
ing that tp^ =< 1 3/(cc) ll 3/(y) >. The superscript is omit- 
ted hereafter. The stationary condition is transformed 
into the following equation, 



d 3 q 



A(g) 



f -{p-qf ^{q)-q 



where C 2 = {N 2 - 1)/(2JV C ) = 4/3 is the quadratic 
Casimir operator for the color SU(3) group. This is noth- 
ing but the Schwinger-Dyson equation for an effective 
quark mass A(p) in the ladder approximation, which is 
the same equation as that obtained by the previous au- 
thors 0- Of course we may take a variation of Delta(p) 
instead of f(p). The stationary condition with respect 
to A(p) leads to 



<fi(p) = 4£ 



A(p) 



&{P)-P 2 



(18) 



Substituting this equation into Eq.(13), we get the same 
Schwinger-Dyson equation. 

If the effective quark mass function is replaced by the 
mean value, we can obtain the partition function at once. 
Moreover from the partition function, the thermodynam- 
ical potential Q is given by 



0(T,/i) = ~S<°> 



T 



d 3 P 

(27T) 5 



log(A 2 (p)-p 2 ) 



+ 



A 2 b) 



A 2 (p) 



P 2 



(19) 



The thermodynamical potential 0(T, /i) will be a key 
function in later discussions on the chiral phase transition 
and the latent heat. 



IV. PHASE DIAGRAM 

We are now in a position to solve the Schwinger-Dyson 
equation (S-D eq.) numerically and determine the phase 
diagram of the quark matter. Let us now adopt the mod- 
ified running coupling of the one-loop order (~g) . It is 
obtained through the study using the operator product 
expansion and renormalization- group equation in QCD. 
The explicit expression is given by 



9 



tip) 



2?r 2 a 



log[(-p 2 



P|)/A 2 



(20) 



QCDl 



where a = 6(72/(11 — 2Nf/3) — 8/9 and the parameter 
Pr is introduced in order to regulate the infra-red diver- 
gence. Here we have assumed that the running coupling 
is independent of chemical potential as done in Ref.jjj. It 
should be noted that the asymptotic freedom in the deep 
Euclidean region is satisfied by this running coupling. 



Instead of solving the S-D equation exactly, we take a 
variational procedure according to Ref. Q : we introduce 
the following trial mass function (A(p)) obtained through 
the study of T = /x = case: 

A W - -^tH(~P 2 +Pk)/A? Q c D ]*~\ (21) 

P ^ Pr 



where a is an order parameter of the chiral symmetry. 
Then the S-D equation with the running coupling is lead 
to 



A(q) = 3C 2 /T 1 E 



d 3 p g 2 (g,p) A(p) 
(2tt) 3 -(q-p) 2 A 2 (p)-p 2 



(22) 



Since — (q—p) 2 in the denominator is a gluon momentum, 
the most natural form of the running coupling would 
be g 2 (— (q — p) 2 )- However, the momentum dependence 
would bring about many difficulties from its angle de- 
pendence in actual numerical calculations. Assuming 
that the effect of the angle dependent part —2qpcos9 in 
— (q — p) 2 is negligible on average, we approximate the 
running coupling [TJ|: 

g 2 (q,p) 



9 2 {-{l 2 +P 2 ))- 



(23) 



Since we use the variational procedure with respect to 
a instead of the self-consistent iteration method, it is 
sufficient to consider the equation with the lowest value 
of -q 2 : —q 2 = ir 2 T 2 . Thus the S-D equation we ought 
to solve is 



1 



^ 2 T 2 + p\ 
32T 



[\og{(n 2 T 2 +p 2 R )/A 2 QCD }\ 
1 



£ / rflPlP 2 ^" 

„ J l0g[(7T 



1 



^ + u J l+^+p 2 R )/K 2 QCD ] 
A{P) (24) 



K - ttT) 2 + p 2 A 2 {p)-p 2 

Since the left hand side is independent of tr, it is rather 
easy to seek the values of a satisfying this equation. 

The parameters of our model are Aqcd = 738 MeV 
and tp = log(p R /A-Q CD ) = 0.1. These values are the 
same as those used in Ref.|7j, where the trial mass func- 
tion is also the same as ours. 

The above equation has three types of solutions in ac- 
cordance with the values of the temperature T and the 
chemical potential /i: 

(1) a = 

(2) a = , ai (0 < a x ) 

(3) a = , CTO , <7l (0 < (T < £71 ) 

In the case (1), a — is always a minimum and the ther- 
modynamical potential is monotone increasing with a. 
In the case (2), a = is a maximum and oi a minimum 
so that Q is decreasing (a < a\) and increasing (a > ui). 
In the last case, a = and a\ are minima and cto is a 
maximum: Q is increasing (cr < cro and a > <7i) and 
decreasing (a < a < a\). Moreover the case (3) is dis- 
criminated between the case of £1(0) > fl(ai) (3~) and 
that of O(0) < fi(cri) ((3+)). 
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FIG. 1: The schematic view of the phase diagram on the T — p 
plane based on our numerical results. The point P denotes 
the tricritical one. 



Now we consider behavior of the thermodynamical po- 
tential when the temperature increases (the chemical po- 
tential is fixed). The numerical result shows three pat- 
terns: 

(a) (2) - (1) or (3~) - (1), (p < p P ) 

(b) (3") - (3+) - (1), (fip < p < po) 

(c) (3+) - (1) or (1). (p <p) 

In the pattern (a) , the order parameter a changes from <j\ 
to at T = T c continuously so that the transition is of the 
second-order. On the other hand in (b), a changes from 
o\ to at T — T c discontinuously. This abrupt jump oc- 
curs under the condition: fi(T c , p, a\) = fl(T c , p, 0). This 
equation determines the critical temperature T c and this 
transition is of first-order. The last pattern does not give 
any phase transitions and corresponds to the supercool- 
ing. The fip m the above equations is a critical point 
between such two transitions and the point {p Pl T P ) is 
called a tricritical point. 

Numerical results are shown in Fig. 1, which is a phase 
diagram for the quark matter. The closed circles mean 
the second-order phase transition and the open ones the 
first-order one. The point "P" stands for the tricriti- 
cal point. We have a second-order phase transition at 
T c = 200 MeV with p = and a first-order phase tran- 
sition at p c = 400 MeV with T = 0. The position of 
the tricritical point "P" is at (T P , p P )=(105 MeV, 300 
MeV). These numerical results are consistent with those 
by the previous authors. 



V. LATENT HEAT 

In this section, we will calculate the latent heat gen- 
erated across the first-order phase transition. For this 
purpose, we need to know the entropy density S and the 
particle number density p first. We note the following 
relation: 



FIG. 2: Temperature dependence of AQ at p = 350 MeV. 
This is a case of the lst-order phase transition. The curve 
passes through the horizontal line at the critical point T c 
and the slope of the tangential line gives the latent heat by 
Eq.(30). 



Accordingly, we can calculate S and p from the thermo- 
dynamical potential as follows, 

S = -8f> P =~W (26) 

Let us now introduce the difference between the thermo- 
dynamical potential in the symmetry broken phase and 
that in the symmetry one, 

AO(T c , p) = fi(T c , M, vi) - n(T c , p, 0). (27) 

Similarly, we define the difference between the entropy 
density in the broken phase and that in symmetry one: 
AS(T c ,p) = S(T c ,p,ai) - S(T c ,p,0). 

We can calculate AS(T C , p) from the partial derivative 
of Ail with respect to the temperature T, 



AS(T C)/ ,) = -(^P) (T C)M ). 



(28) 



Now the latent heat can be derived from this entropy 
difference as follows: 



Q l = -T c AS(T c ,p). 



(29) 



For example, we have drawn the Af2(T, p) as a function 
of the temperature with p = 0.35 GeV in Fig. 2. Then 
we can evaluate the latent heat from the slope of the 
curve at the critical point AQ(T c ,p) = 0. The result 
of the latent heat per particle thus obtained is shown in 
Fig. 3 as a function of the chemical potential p. It tells 
us that the latent heat is not given off in the second- 
order phase transition (p < p P = 0.30 GeV). The latent 
heat emerges abruptly with large value (<~ 250 MeV) 
near the tricritical point (pp = 0.30 GeV). When the 
chemical potential is increased, it decreases and vanishes 
at p = 0.40 GeV. In order to realize this behavior, let us 
consider the following relation, 



dn = -SdT~ pdp. 



(25) 



(30) 
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FIG. 3: The latent heat Qi as a function of the chemical 
potential fi. 

which is derived from Eqs.(28) and (29). When ijl is in- 
creased in the region of the first order transition, the 
figure of Afl(T) becomes shallow as a function of T (see 
Fig.2): This means that the derivative dAfl(T c )/dT c de- 
creases with increasing fi. Therefore the latent heat de- 
creases with increasing [i. Our numerical result shows 
that the latent heat has a singularity near the tricritical 
point and plays an important role to identify the tricrit- 
ical point. 

VI. SUMMARY AND DISCUSSIONS 

We have studied the latent heat in the chiral phase 
transition in this paper. We have taken the QCD-like 
gauge field theory as an effective theory and have per- 



formed the numerical calculations using the mean field 
(WKB) approximation. We have calculated the latent 
heat generated in the first-order phase transition. The 
latent heat is enhanced near the tricritical point and de- 
creases as the chemical potential increases. However, it 
should be pointed out that more accurate numerical sim- 
ulations will be necessary for the detailed behavior of the 
latent heat near the tricritical point. 

We must now note the existence of the color super- 
conducting state which has been disregarded in our nu- 
merical calculations [23]-[24|. The color superconducting 
state is considered to be able to exist in the region at low 
temperature and high density apart from the tricritical 
point. Hence, we can not apply our results directly on 
the actual physical phenomenon, in which color super- 
conductivity is involved. However, the phase transition 
on the color superconductivity is of the second-order so 
that it does not give off the latent heat originated on this 
transition. 

It is pointed out that the values of T and /i realized 
in high energy heavy-ion collisions may be close to the 
tricritical point 0I"E3- Therefore, it may be possible 
to observe some signals originated in the generation of 
the latent heat because of the fact that the latent heat 
becomes large near the tricritical point. 
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